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ABSTRACT 

The topic of this article is weak cosmic shear tomography where the line of sight-weighting is 
carried out with a set of specifically constructed orthogonal polynomials, dubbed TaRDiS (To- 
mography with orthogonAl Radial Distance polynomial Systems). We investigate the proper- 
ties of these polynomials and employ weak convergence spectra, which have been obtained by 
weighting with these polynomials, for the estimation of cosmological parameters. We quan- 
tify their power in constraining parameters in a Fisher-matrix technique and demonstrate how 
each polynomial projects out statistically independent information, and how the combina- 
tion of multiple polynomials lifts degeneracies. The assumption of a reference cosmology is 
needed for the construction of the polynomials, and as a last point we investigate how errors in 
the construction with a wrong cosmological model propagate to misestimates in cosmological 
parameters. TaRDiS performs on a similar level as traditional tomographic methods and some 
key features of tomography are made easier to understand. 
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1 INTRODUCTION 

Weak gravitational lensing by the cosmic la r ge-scale structure 
JBlandford et alj Il99ll; j Schneider et al.l 1 19921 : Iseitz eTal] 1 1994 
ISeitz & SchneideJ [19941 : Kamionko wski et alj Il998l) is regarded 
as a very promising way for investigating the properties of the 
cosmic matter distribution and the measurement of cosmological 
parameters. The prim ary tool is the angular spectra of the weak 
lensing c onvergence i Jain & SeliakI 119971 : IHu & TegmarM 1 19991 : 
IHu& Wh ite 2001; Hu & JairJboCT . which has, in the context of 
dark energy cosmologies with adiabatic fluctuations in the dark 
matter distribution, the potential to deliver parameter constraints 
competitive with those from the cosmic microwave background. 

The weak cosmic shearing effect of the large-scale structure 
has b een detected by four independent research groups a decade 



ago (Van Waerbekeetal.! 200i 



le pendent research groups a decade 
03: iKaiser et alj l2000l : iBacon et al.1 



I2OOO; : .Wittman et aL .20001) and now parameter constraints de 
rived from weak lensing spectra coincide well with those from the 
CMB, with a small tension concerning the parameters Q.„, and erg. 
(Kilbinger et al. 2009, 2010). These topics are reviewed i n detail by 
fBarte lmann & Schneideil ( 1200 ih and lBartelmaiii] i2010J lbl). 

Together with the property of weak lensing to map out fluc- 
tuations in the cosmic matter distribution in a linear way, tomo- 
graphic methods have been shown to offer superior precision in 
the determination of cosmological parameters (Hu 1999, 2002a. b; 
iTakada & White! |2004 iHannestad et alj|2006l). in particular con- 



cerning the p r operties of dark energ y dHuterer & Tumei 
iHutereJlIooa : iHollenstein et alj 120091 : lHutereJl2010l : iLee 
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lAmara & Kitchinal201lh . The idea of tomography is a division of 
the galaxy sample used for shear estimation into a number of bins 
in distance. This splitting is the reason for two important advan- 
tages of tomography: The weak lensing convergence is a line of 
sight-integrated quantity of which one measures the angular fluc- 
tuation statistics. Angular convergence spectra, however, convey 
less information compared to the fluctuations statistics of the three- 
dimensional matter distribution because in the projection process, 
a mixing of scales is taking place. Secondly, cosmological param- 
eters can have different influences on the amplitude of the weak 
lensing convergence at different redshifts, which would be aver- 
aged out. Combining convergence spectra measured for each sub- 
set of galaxies are able to alleviate these these two problems. Re- 
lated developments include shear ratio measurements, where the 
same tidal fi elds are measured with galaxy populat ions at differ- 
ent redshifts Uain & Tavlonl2003l : iTavlor et alJbOOTT) , cross corre- 
lation cosmography, where part of the galaxies are used for con- 
structing a template on which one observes th e shearing effect of 
more distant galaxies ( IBernstein & Jainll2004h . and finally three- 
dimensional weak shear methods, which provide a direct recon- 
struction of the fluctuations of the cosmic densit y field dCastro et alj 
'2005':'ffeav ens et al.l2 00?:'Kitching e t alj20Tlh . Common to these 
methods are nonzero covariances between measurements of differ- 
ent redshifts. This is a natural consequence that lensing is caused 
by the large-scale structure between the lensed objects and the ob- 
server, so that the light from different galaxy samples has to trans- 
verse partially the same large-scale structure for reaching the ob- 
server 

Our motivation is to revisit weak lensing tomography and to 
construct a weighing scheme for the galaxies as a function on dis- 
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tance, such that the covariances between spectra estimated from 
differently weighted galaxy samples has a particular simple, diago- 
nal shape. For this goal, we construct a set of orthogonal polynomi- 
als in distance, which provides a diagonalisation of the covariance, 
and which, when employed in the derivation of the convergence 
spectrum, projects out statistically independent information about 
the large-scale structure. 

After a brief summary of key formulae for cosmology, struc- 
ture formation and weak lensing in Sect. |2l we describe the con- 
struction of TaRDiS-polynomials and investigate their properties 
in Sect. [3] Their statistical and systematical errors are evaluated 
in Sects. |4] and [5] respectively, and we summarise our results in 
Sect. |6l The reference cosmological model used is a spatially flat 
vi'CDM cosmology with Gaussian adiabatic initial perturbations in 
the cold dark matter density field. The specific parameter choices 
are Sl^ = 0.25, erg = 0.8, Ho = 100 h km/s/Mpc, with h = 0.72, 
nj = 1 and flj = 0.04. The dark energy equation of state is set to 
w = -0.9 and the sound speed is equal to the speed of light (c, = c) 
such that there is no clustering in the dark energy fluid. 



2 COSMOLOGY AND WEAK LENSING 

2.1 Dark energy cosmologies 

In spatially flat dark energy cosmologies with the matter density 
parameter Q,,,, the Hubble function H{a) = dlna/dt is given by 



(1) 



H^(a) Q.„, I f' 

— 5- = ^- + (l-n„)exp 3 dlna[l+vKa)] 

with the dark energy equat ion of state w(a), for which we use the 
common parameterisation JChevallier & Polarskil200ll) 



w(a) = Wq + (1 - a)Wa, 



(2) 



where wo = - 1 and w^ = would correspond to ACDM. Comov- 
ing distance x and scale factor a are related by 



X = c 



r' da 

1 a^Hia)' 



(3) 



such that the comoving distance is given in units of the Hubble 
distance ;t'/; = c/Hq. 



2.2 CDM power spectrum 

The linear CDM density power spectrum P(k) describes statisti- 
cally homogeneous Gaussian fluctuations of the density field S, 



{6(k)6(k'r) = {InfSoik - k')P{k). 



(4) 



It is composed from a scale invariant term oc k'' and the transfer 
function T{k), 



P(k)cck"'T^(k). 



(5) 



In low-Hm cosmologies T(k) is approximated with the fit proposed 
bv lBardeenetaLHl986l) . 

ln(l -l-2.34<7) / , , ,^-i 

T(q) = ' ^ 3^ (1 + 3.89^ + (l6Aqf + (5A6qf + (GJlq)") ' 

where the wave vector k = qT is given in units of the shape param- 
eter F. F determines the peak shape of the CDM spectrum and is 
to first order given by F = il„,h, wi th small corrections due to the 
baryon density Hj, l lSugivamalll995f) . 



F = Q„,h exp -Qi 1 -1- 



V2^ ]] 



(7) 



The spectrum P(k) is normalised to the variance os on the scale 
R = 8 Mpc/h, 



^2 ^ fdlnk ,,3j,,,^„„2, 
" J 2; 



In^ 



eP(k)Wl{k) 



(8) 



with a Fourier transformed spherical top hat filter function, WR{k) = 
3ji (kR)/(k R). /V(x) is the sp herical Bessel function of the first kind 
of order { JAbramowitz & Stegua. 19721 



2.3 Structure growth with clustering dark energy 

Linear homogeneous growth of the density field, S(x,a) = 
D+(a)d{x,a = 1), is described by the g rowth function D+(a) , 
which is the solution to the growth equation 1 Turner & White! 19971 : 
I Wang &Steinhardtll 19981: iLinder & Jenkin3l2003h . 

^DJa) + - (3 + ^^1 ^D^a) = -^fi,„(a)£>+(a). (9) 

da-^ a\ dlna / da 2a^ 



Nonlinear structure formation enhances the CDM-spectrum P(k) 
on small scal es by a f actor of == 40, which is described by the fit 
suggested bv lSmith et al. (2003) which is gauged to n-body simu- 
lations of cosmic structure formation. 



2.4 Weak gravitational lensing 

The weak lensing convergence k provides a weighted line- 
of- sight average of the matter density 6 (fo r reviews, see 
iBartelmann & Schneide3l200ll : lBartelmanrJl2010bh 



Jo 



dX W,(x)S, 



(10) 



r dx'q(z)p;^^.(n) 

Jy 4^' X' 



with the weak lensing efficiency Wi^ix) as the weighting function, 

w.Or) = J^— GOr)^, withGCt) 

qiz) denotes the redshift distribution of the lensed background 
galaxies, 

I /n I .1 I 

(12) 



q{z) = ?o I — I exp 



d, with 1 = 5If- , 



and is the approximate forcasted distribution of the EUCLID mis- 
sion ( Amara & Refregier 2007). zo has been chosen to be ^ 0.64 
such that the median of the redshift distribution is 0.9. With these 
definitions, one can carry out a Limber-projection ( lLimbenll954h 
of the weak lensing convergence for obtaining the angular conver- 
gence spectrum Cuit), 



CM 



% Wl{x)P{k = tlx\ 



(13) 



lo X 
which describes the fluctuation statistics of the convergence field. 



3 TOMOGRAPHY WITH ORTHOGONAL 
, (6) POLYNOMIALS 

3.1 motivation for tomography 

Weak cosmic shear provides a measurement of the weighted, line of 
sight integrated tidal shears althoug it is more convenient to work in 
terms of the weak lensing convergence and the cosmic density field. 
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which have statistically equivalent properties. Weak lensing there- 
fore provides an integrated measurement of the evolution of the 
cosmic density field weighted with the lensing efficiency function. 
Being a line of sight-averaged quantity, the weak lensing conver- 
gence is statistically not as constraining as the full 3-dimensional 
density field, which is caused by the mixing of spatial scales in the 
Limber-projection. Additionally, parameters which enter the model 
in a nonlinear way can cause different effects on the weak lens- 
ing signal at different redshifts. This sensitivity would be averaged 
out in a long-baseline line of sight integration and the information 
would be lost. 

The solution to this problem are tomographic methods: split- 
ting up the lensing signal from different distances allows a larger 
signal strength and a higher sensitivity with respect to cosmologi- 
cal parameters because the line of sight-averaging effect is reduced. 
This comes at the cost of a more complicated covariance between 
the spectra measured, and a higher shape noise for each of the spec- 
tra, because the effective number of galaxies is reduced. There are 
different methods for exploiting evolution of the cosmic density 
field and the lensing sensitivity along the line of sight. Since the in- 
ception of tomography for investigating weak lensing convergence 
spectra (Hu 1999, 2002b) and bispectra (Takada&Jain 2004), 
there have been many de velopments leading u ltimately to cross- 
correlation cosmography JBemstein & Jainll 20041 , where a lensing 
template is derived from data and used to predict the weak shear 
signal on background sources, and to 3-dimensional weak shear 
methods, which are an unbinned, direct mapping of the cosmolog- 
ical d ensity field JCastro et al.ll2005l : lHeavensll2003l : [Kitching et al] 

lioTIh . 

Common to these methods are strong covariances between 
spectra, which we aim to avoid by employing a set of specif- 
ically designed polynomials, which perform a weighting of the 
observed ellipticities in their distance, in a way that convergence 
cross-spectra computed with two different polynomials vanishes. 
This ensures that every polynomial projects out information about 
the deflection field which is statistically independent, and gener- 
ates a covariance matrix of a particular simple shape. In a sense, 
TaRDiS is using a non-local binning of the ellipticities motivated 
by the anticipated weak lensing signal, which corresponds to the 
non-local nature of weak lensing. 



with the weighted convergence spectrum Sjjit) for homogeneous 
and isotropic fluctuations of /c;. 



s.M) - 



"" ^w,(x)WJ(x)P(k = e/x)■ 



im 



Jo X 

This expression provides the definition of a scalar product between 
the weighting functions p,(x) and Pjix) contained in Wiix) and 
y^ j(x), respectively. 



<p„p,> = 5„(a 



(19) 



which exhibits the necessary properties of being positive defi- 
nite {(pi,Pi) > and (Pi,Pi) = Q <r^ Pi = 0), symmetric 
i(PhPj) = (Pj, Pi)) and linear in both arguments. Starting point 
for the construction of orthogonal polynomials which would diag- 
onalise S ij{{), S ij(C) oc Sij, is the family monomials 



pW 



— ■ 



(20) 



where ;i^„odc can be chosen for placing the node of the first polyno- 
mial in a sensible way, in our case it is set to the median value of 
the galaxy redshift distribution, converted into comoving distance. 
The monomials p'j(x) are subjected to a Gram-Schmidt orthogonal- 
isation procedure, with the initial condition 



Pq(x) = p'q(x) = 1 
and, iteratively, 

Pi(x) = p'i(x) - 2 



(P'rPj) 
(Pj'Pj) 



Pj(x)- 



(21) 



(22) 



It should be emphasised that the orthogonalisation procedure needs 
to be repeated for every multipole £, and that we have omitted an 
additional index i of the polynomials p,Of ) for clarity. For illustra- 
tion the polynomials are normalised using the norm induced by the 
scalar product defined in eqn. l |19t . 



Piix) ■ 



Piix) 



(23) 



which appears in the denominator of eqn. j22t in a natural way. 
Clearly, the scalar product 5,^(0 is equal to the convergence spec- 
trum Ck({) for i = j = 0, in which case po = 1- 



3.2 construction of orthogonal sets of polynomials 

We introduce a weighting of the distance distribution of the galax- 
ies q(j(;) = q(z)dz/dx = q(z)H(z) with a function Piix), such that 
the lensing efficiency function reads 



3fim D+ 

^xi " 

with 



(14) 



f 



G,{x)= I Ax'q(x)p,ix)^^^- (15) 

•Jx X' 

A line of sight-weighting of the weak lensing convergence with a 
function PiOf) yields the weighted convergence field /c,. 



Jo 



Ax iv,Cr)5, 



(16) 



from which one obtains for the variance between two Fourier- 
modes Ki{E), 



{k,{{)k]{{')) = (2nfSo(f-nS,-j(f), 



(17) 



3.3 properties of orthogonal polynomials 

Fig. [T] shows the polynomials Piix) at a fixed multipole order of 
{ = 1000. They exhibit a sequence of oscillations roughly at the 
positions where the previous polynomial assume their maximal 
values. When computing the polynomials for the nonlinear CDM 
spectrum instead for the linear one, the sequence of oscillations 
is shifted to smaller comoving distances. The reason for this shift 
is the larger amplitude of the nonlinear spectrum P(k) for small 
k — Cjchi, which causes a larger contribution of the weak lensing 
spectrum to be generated at high x- 

The variation of the polynomials p,Of) with distance x and 
multipole order I is depicted in Fig. |2l in this case for the poly- 
nomial ( = 9. Again, small variations with multipole order E are 
present, and the oscillations when varying x are clearly seen. 

The orthonormality relation {phPj) of of the polynomials 
Piix) for the lensing signal as well as for the galaxy shape noise is 
given in Fig. [3] The polynomials are constructed to be orthogonal 
by the Gram-Schmidt procedure, but numerical noise is collected 
in the iterative process, such that the orthogonality relation is better 
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Figure 1. Orthogonal polynomials Pi(x), i = ... 9, as a function of comov- 
ing distance x constructed with the Gram-Schmidt algorithm for the weak 
tensing spectrum at f = 10^, with the lowest order polynomials on top, for 
both a linear (solid line) and a nonlinear (dashed line) CDM spectrum P(k). 
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Figure 2. Orthonormal polynomials Pi(x) as a function of multipole order 
i and comoving distance X- The order of the polynomial has been fixed to 
i = 9 and was constructed for a nonlinear CDM spectrum. One can see the 
slow variation of the Pi(x) polynomials with multipole order C. 



fulfilled at small ; compared to larger i. Deviations from {/?,, pj) = 
for ; t j are of the order of ~ 10 '"^ for low-order polynomials, but 
increase to values of ~ 10"'* at high order. This deterioration in or- 
thogonality is a known drawback of the Gram-Schmidt procedure, 
in particular when dealing with sets of functions instead sets of vec- 
tors, which means that there is larger numerical noise in the evalua- 
tion of the scalar products, but is acceptable for the number of basis 
polynomials we are using in this work. One can already notice the 
main difference between the different approaches in tomography: 
Whereas the noise covariance would be diagonal in classic tomog- 
raphy with a very complicated structure of the signal covariance, in 
our case the shapes of the signal and noise covariance matrix are 
interchanged. 

Fig. |4] gives an impression of the lensing efficiency function 
Wj(x) modified by the polynomials Pj(x) on an angular scale of 
£ = 1000, in comparison to that of weak shear without tomography, 
Wo(;i') = Wit(j() for i = 0. It is quite interesting to see how the 




Figure 3. Orthonormality relation for the polynomials Pi(x), at f = 10^ in 
logai'ithmic representation, for the signal covariance 5 ,j (left panel) and the 
noise covariance Njj (right panel), for the EUCLID survey characteristics, 
in units of the shape noise <T-/n. 




comoving distance ;j^ in Mpc/ft 

Figure 4. Lensing efficiency functions Wi(x), i = ... 9, as a function of 
comoving distance^, for both linear (solid line) and nonlinear (dashed line) 
CDM spectra. The lensing efficiency functions have been constructed for 

e= 10^. 



seemingly messy functions Wjix) disentangle and approach zero at 
large distances. 

Finally, polynomial-weighted weak lensing spectra 5„(f) = 
{pi,Pj) are shown in Fig.|5] The spectra drop in amplitude, which 
is mostly an effect of the absence of normalisation, and there are in 
fact differences in shape, when higher-order polynomials are used 
(compare Fig. lAll in the appendix, where all spectra are scaled such 
that they assume the same value at a certain multipole.). By con- 
struction, these spectra provide statistically independent informa- 
tion of the cosmic large-scale structure. In the next sections, we 
will investigate statistical bounds and systematical errors on cos- 
mological parameters, when the information from different spectra 
is combined. 



4 STATISTICAL ERRORS 

In this section we show how we construct covariance matrices for 
the polynomial-weighted spectra, compute the Fisher-matrix and 
derive statistical errors on cosmological parameters, and investigate 
the signal strength of the weak shear signal, all as a function of the 
number of polynomials used. 
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Figure 5. Weak lensing spectra Salt), i = . . .9, weighted with orthogo- 
nal polynomials piix), as a function of multipole order £, for both a linear 
(solid line) and a nonlinear (dashed line) CDM-spectrum P(k). The spectra 
decrease in amplitude with increasing polynome order ;' for unnormalised 
polynomials. 



4.1 variances of weighted ellipticities 

For deriving the expressions for the signal covariance and the noise 
covariance needed in forecasting statistical and systematical er- 
rors, we derive expressions for the mean and the variance for a 
weighted set of ellipticities. This derivation is done for a discrete 
set of weighting coefficients w„, and then generalised to the con- 
tinuous case, where the weighting is done with a polynomial p,(^). 
The distribution p(e)de of ellipticities e is assumed to be Gaus- 
sian, with a zero mean, variance a-^ and the ellipticities are taken 
to be intrinsically uncorrelated (i.e. intrinsic alignment effects are 
discarded, see Schafer 2009|, for a review). 

From a measurement of a set of ellipticities 6,„ drawn from the 
parent distribution p(e)de on can estimate the shear by computing 
the weighted mean e, and the expectation value of the weighted 
mean (e) if the drawing of ellipticities is repeated, 



<e>: 



^0, 



(24) 



2jm ^^iii 2-tm ^''" 

which vanishes if the mean ellipticity vanishes, {e„,) = 0. For the 
variance {€'} one obtains under the assumption of intrinsicalUy un- 
correlated ellipticities ,(e,„e„) = cr^S„,„, 



C^) 



1 



■ /^W„W„{(E„,£„) 



■Z' 



(25) 



which reduces to the classic Poissonian result if w,„ is either or 1 



(tr) 



— ^ with 

A' 



A': 



Z' 



(26) 



because w^, = w„, and A' is defined as the effective number of 
ellipticities in the sample. The cross-variance for two dilTerent sets 
of weights w„, and v„ is given by: 



Ce') 



■z 



w,„v„, 



Urn '^m Z;i ''" 

In the continuum limit we make the transition 
^ . . . ^ n I dxq(x) ■ ■ ■ 



(27) 



(28) 



with the unit normalised galaxy distance distribution q(x)dx- The 
discrete weights w,,, and v„, will be replaced by the set of poly- 
nomials Piix) and Pjix)- For conserving the normalisation of the 
weighted galaxy distance distribution q{x)Ax^ we normalise the 
polynomials 



P.(X) 



P.(X) 



j dx q(x)Pi(x)' 



such that the weak shear spectrum 5, y(^) becomes 



Sij = (PhPj) 



" % W,{x)Wj(x) P(k = dx). 



which is diagonal by construction, and the noise covariance. 



N.i 






dx q(x)Pi(x)Pj(x) 



(29) 



(30) 



(31) 



with the shape noise cr^ and the mean density of galaxies h per 
steradian, for which we substitute the numbers 0-,: = 0.3 and 
h = 40/arcmin^ projected for EUCLID. It can already be seen 
from the expression for Njj that the weighting of data with high- 
order polynomials piix) will be noisy: qix) is a slowly varying 
function and the integrals jdx qix)Pi(x) in the denominator will 
assume small values if the polynomial is rapidly oscillating. This 
will ultimately limit the order of the usable polynomials. Again, 
for i = j = 0, the standard Poissonian expression A'oo = crl/h is re- 
covered due to the normalisation of q(z), as well as the weak shear 
spectrum 5oo(^) = Cn({)- 



4.2 Fisher-analysis 

The likelihood function for observing Gaussian-distributed modes 
Ki(t) of the p,(^)-weighted weak lensin g convergence for a 
given parameter s et x^ is defined as (see iTegmark et alj 119971 : 
ICarronetalj201ll) : 



JK^mx.) 



V(2;r)'Vdet(C) 



exp(-L,(e)Ci/(e,e')K'j(r) 



(32) 



with the covariance Cij(t,t') = {Ki(t)K'(t')) which is diagonal 
in £ for homogeneous random fields. The ;^'^-functional, X 0^ 
exp(-;^'^/2), can be obtained from the logarithmic likelihood and 
reads: 



X' 



2 = ^ tr [in C-^C^'o], 



(33) 



with the definition of the data matrix D^j = Ki(l)Kj{t), using the 
relation In det(C) = tr ln(C) and discarding irrelevant multiplicative 
prefactors. The second derivatives of the ^^ -functional with respect 
to cosmological parameters x^ evaluated at the point of maximum 
likelihood yields the Fisher matrix 



d' X- 
dxudxy 2 



V-i 2f +1 
> ^ — tr 
V 2 



d d 

lnC--lnC 

OXn OXy 



(34) 



with multiplicity 2t + I because on each angular scale t there are 
2^+1 statistically independent ra-modes. The covariance Ctj = 
S ij + Njj can be split up into the signal covariance S,j, 



C Soo(0 = CM 



s,m 







^ qq\^) , 



(35) 



and the noise covariance Nn 
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Figure 6. Derivatives -Jtr (da„, In C) (solid lines) and -Jtr (5„ In C)" 
(dashed lines) as a function of multipole order f and cumulative polynomial 
order q. The derivatives have been weighted with the EUCLID covariance. 



N.AO 



( iVooW = o-l/n 



No,({) 



N^(() 



N,,({) ) 



(36) 



We will work in the limit dS ij/dxij » dNfj/dxij which is well justi- 
fied in our case. 

Fig. [6] gives an example of how the combination of multiple 
line of sight-weighted measurements helps to avoid weaknesses in 
the sensitivity towards cosmological parameters. We plot the sensi- 
tivity 



ainc 

dxu 



dF„. 



2£+\ de 



(37) 



i.e. the ratio between the derivative of the spectrum with respect to 
a cosmological parameter and the covariance of the measurement, 
which corresponds to the contributions to the diagonal elements of 
the Fisher-matrix per f-mode. The sensitivites of the measurement 
with respect to the cosmological parameters Q„, and w are shown, 
which exhibit singularities at certain multipoles for a classical non- 
tomographic measurement. This happens when a parameter affects 
a certain /"-range with a different sign than others and the spec- 
trum pivots around this value when varying a parameter. Clearly, 
angular scales in the vicinity of that pivot scale do not add much 
sensitivity to the Fisher-matrix. The inclusion of a second polyno- 
mial, however, avoids this: For ^ > 1 on, singularities are lifted 
and the derivatives assume larger values with increasing numbers 
of polynomials, although this effect saturates for very large num- 
bers of polynomials. At large multipoles, the influence of the shape 
noise can be observed, which causes the derivatives to drop rapidly, 
and do not add significant sensitivity to the Fisher matrix on scales 
larger than £ ^ 3000. 



4.3 statistical errors 

From the Fisher-matrix F^v one can obtain the Cramer-Rao errors. 



(p-%^ 



(38) 



and the 2-dimensional marginalised logarithmic likelihood Xm 
around the fiducial model x'. 
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Figure 7. Statistical errors on the estimates of the cosmological parameters 
Cl,„ (circles), erg (squares), h (lozenges), (i., (triangles, pointing up) and w 
(triangles, pointing down) resulting as Cramer-Rao errors from computing 
the Fisher-matrix for the weak lensing spectrum SaiC), as a function of 
cumulative polynomial order q, for EUCLID survey characteristics. The 
maximum multipole consideres is f^ax = 1000 (dotted lines) and /"i^ax = 
3000 (solid hues). 
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Fig. |7] summarises statistical errors on cosmological parame- 
ters resulting from the Fisher-matrix analysis, cumulative in q and 
for £ = 1000 as well as for £ = 3000. As expected, statistical errors 
drop with larger number of polynomials and are smaller if more 
multipoles are considered. A very fascinating feature of the plot is 
the approximate scaling of the error with the inverse root of the 
number of polynomials, cr oc 1/ -yjq, which one expects from Pois- 
sonian arguments because each spectrum 5„(f) adds statistically 
independent information. 

2-dimensional marginalised likelihoods for all parameter pairs 
are shown in Fig. [8] for £ = 1000 and £ = 3000 as the maximum 
multipole considered. In the confidence contours, up to 10 poly- 
nomials were combined. Clearly, the uncertainty in cosmological 
parameters decreases with larger numbers of polynomials used, as 
well as for an increased multipole range. Additionally, the param- 
eter degeneracies change a little when more polynomials are used, 
which is caused by differing sensitivities of the lensing signal with 
distance. A 3-dimensional view of the marginalised likelihood of 
the parameters fi„,, erg and w is given in Fig.|9] It illustrates nicely 
the nested Icr-ellipsoids which become smaller for larger number 
of polynomials combined in the measurement. Comparing the per- 
formance of the TaRDiS-polynomials to other tomographic meth- 
ods shows that they operate on very similar levels of performance, 
perhaps with a small advantage for the cosmological parameters Q.,„ 
and w. This is due to the mechanism that the values in the Fisher- 
matrix are maximised if the covariance becomes diagonal due to 
a match between the true cosmology and the assumed cosmology 
used for constructing the polynomials. This should be taken with 
a grain of salt, however, because part of the statistical error would 
be transported to the systematical error budget. The impact of start- 
ing off with a wrong cosmological model for construction of the 
TaRDiS-polynomials on the estimation of parameters will be the 
topic of the next chapter. 



tomography with polynomials 7 



0.26 



a 



0.25 



0.24 




0.78 



0.8 
OS 



b 0-8 



0.78 








0.65 



0.7 

h 



0.7 



-5; 0.7 



0.65 





0.95 



1 



1.0 



0.95 



-1 



-0.9 
W 



-0.8 



Figure 8. Consfraints on the cosmological parameters CI,,,, erg, h, «, and w from EUCLID using tomography with orthogonal polynomials. The ellipses mark 
Icr confidence regions and decrease in size with increasing cumulative polynomial order (^ = in blue to ^ = 9 in green), for ^,nax = 1000 (dotted lines) and 
f,^^x = 3000 (solid lines). 



4.4 signal to noise-ratio 

Analogous to the definition of the Fisher-matrix we construct the 
cumulative signal to noise-ratio S, 



c ^ c 



2[+ 1 



tr(c-'5f 



(40) 



from the covariance matrices with the multiplicity 2{+ 1 . The cumu- 
lative signal to noise-ratio S and the differential contribution dS/df 
of each multipole is summarised in Fig. [TO] Clearly, increasing q 
or £ increases the signal up to multipole orders of a few thousand. 
The growth of dE(€)/df oc ylT+l is driven by the reduced cosmic 
variance until the shape noise limits the measurability of the spec- 
tra. Consequently, the integrated signal to noise ratio settles off at 
a few hundred, and adding statistically independent information by 
using new polynomials increases the signal strength by almost an 
order of magnitude and reaches values comparable to the primary 
CMB temperature anisotropy spectrum (Hu 2002b ). These nu m- 
bers correspond well to those derived by Takada & JainI ( 120090 if 
one works in the approximation of a Gaussian covariance - non- 
Gaussian contributions can significantly lower the signal strength. 
This can be expected from eqn. ([40j, as the signal to noise ratio is 
invariant under orthogonal transformations diagonalising either 5,; 
orNij. 



5 SYSTEMATICAL ERRORS 

Naturally, the set of polynomials used for analysing the data needs 
to be constructed for specific cosmology, so the question arises if 
an imprecise prior knowledge of the cosmological model has an 
impact on the parameter estimates from the weighted weak tensing 
spectra. Any incompleteness or imperfection in the model used for 
interpreting the data is going to shift the estimated parameter values 
away from their true values and introduces parameter estimation 
biases. Specifically, we consider three cases: 

Firstly, a time- varying equation of state (w = -0.8-1-0.2(1 -a)) 
when in reality the equation of state of dark energy is constant 
(w = -0.9), i.e. a systematic which is not a degree of freedom of 
the model, and secondly wrongly assumed Q.,,,- and erg-values (0.3 
instead of 0.25 and 0.85 instead of 0.8, respectively) as a strong sys- 
tematic. In addition, we compute biases in the estimation of cosmo- 
logical parameters if the redshift distribution of galaxies substituted 
in the polynomial construction is not the true one. 



5.1 parameter estimation bias 

For the Gaussian likelihood function X o^ exp(-^^/2) for a 
parabolic ;^f^ -functional one can identify x^ - 2f tr In C + C"'£) 
with the covariance Cij and the data matrix £),y. Now, the fit of 
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Figure 9. Icr-ellipsoids in the space spanned by (n,„,cr8,n'), cumulative 
in polynomial order q, from q = 2 (largest ellipsoid) up to g = 9 (smallest 
ellipsoid) with the maximum multipole order set to f = 3000. Observational 
characteristics correspond to those of the EUCLID mission. 
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Figure 10. Signal to noise-ratio S(^) (solid lines) and the contribution dZ/df 
by each multipole (dashed lines) as a function of inverse angular scale C and 
cumulative polynomial order q. The shape noise corresponds to that of the 
projected EUCLID performance. 



the true model C, to the data would give rise to the correct x^- 
functional. 



Xl 



2 tr [in C, + C,-'d], 



whereas the assumption of a wrong model yields 



x] 



2tr[lnC/ + C7'fl], 



(41) 



(42) 



which in general will assume its global minimum at parameter val- 
ues different than;^^^. The distance 6 between the best fit values x, 
of the true model and Xf of the false model will then be the param- 
eter estimation bias. This parameter estimation bias can be quanti- 
fied with a second-order Taylor expansion of the ;^'^-functional for 
the wrong model around the best-fit point x, of the true model (see 
ICabre et al.ll2007l) . 



4(.,) =x}(x,) + 2 ^4(-')^''4 Z ^^/(-')'^''^-^43) 

where the parameter estimation bias vector 8 = Xf - x,. The 
best-fit position Xf of ;^'^. can be recovered by extremisation of the 
ensemble-averaged Ct^>, yielding 



d 



£-; =z- 



d^ 



dx^,dxy 



■Xj Sy, 



which is a linear system of equations of the form 
2_j G^ySy = a,, -» 6,, = y ^(G'')ijyay, 



(44) 



(45) 



where the two quantities G^y and a^, follow from the derivatives of 
the averaged ^^-functional, evaluated at x,, with {D) = C, and the 
multiplicity 2^ -n 1 added for each f-mode: 



Y-i 2^+ 1 



d 



— lnC,(id-C/C,) 



(46) 



This vector reduces to a^ = if C, = Cf (id being the identity 
matrix). Furthermore, 



G„„ 



V-i 2^+1 

= ) tr 

V 2 



^lnC,^lnQ(2C/C,-id) 



(47) 



which simplifies to G^,, = F^y in the case of choosing the correct 
model, such that the parameter estimation bias vanishes. The same 
happens for q = 0, i.e. if only the polynomial poOf ) = 1 is used. In 
the derivation outlined above, the identities 



d 



d 



d 



-— In C = C-'-—C and -— C^' = -C^M -— C C^' (48) 



dx. 



dx, 



dx. 



d 



dx^ 



were used. This formalism is a generalisation of the case of di- 
agonal covariance matric es i Cabre et al.ll2007l . lAmara & Refregien 



l2()08l ; lTaburet et al.l2009l : lMarch et al.ll201 ih . for which it has been 
shown to work well by comparison with results from Monte-Carlo 
Markov chains JTaburet et al.ll201 0). We need to employ this more 
general formalism because the signal covariance 5,^(0 ceases to be 
diagonal if the cosmology used for constructing the polynomials 
differs from the true cosmology. Other examples of non-diagonal 
covariances include 3-dimensional cosmic shear JHeavensluOOSl : 
Kitch ing et alJuOOSh . It should be emphasised that we are only in- 
vestigating the impact of systematics or a wrongly chosen initial 
cosmology on the construction of polynomials and that system- 
atics control plays a very important role in parameter esti mation 

from weak lensing (among others, see 

iHuterer et alj|2006k [Bridle & KindlioOTTlSemboloni et all 



which we are not touching here. 



5.2 systematical errors 

We constructed TaRDiS-polynomials for three wrong assumptions: 
As the first example, we consider a rather small change in the dark 
energy model, namely a time varying gequation of state instead of 
a constant one. These two cosmologies have the same average dark 
energy equation of state, but the degree of freedom of a time vary- 
ing w is contained in the true cosmology. For this case. Fig. fTTI 
shows the estimation bias in the true dark energy model, when the 
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model used for constructing the polynomials had been a time vary- 
ing equation of state with the same average eos-parameter. The fig- 
ure illustrates that such a mistake has a minor impact on the estima- 
tion of parameters, as biases are smaller compared to the statistical 
precision by at least an order of magnitude. Most estimation bias 
can be found in the normalisation erg. 

Biases in the estimation of cosmological parameters if for the 
construction of the polynomials the wrong values for Cl„, and erg 
have been used, are summarised in Fig. [12] This case is a much 
stronger systematic compared to the previous case. There are strong 
biases in particular in n„, and erg, and to a lesser extend in n,, 
whereas h and w are not strongly affected. Furthermore, the biases 
in n,„ and erg are in a direction almost orthogonal to the orientation 
of the degeneracy, indicating a patricularly strong impact. These es- 
timation biases, however, can be reduced by including a larger num- 
ber of polynomials. In that way, a reduction of the bias to values 
similar to the statistical error is possible. 

Finally, biases due to using a convolved redshift distribution in 
the polynomial construction whereas the signal is in reality gener- 
ated by the unconvolved redshift distribution have been computed 
in Fig. [13] As a very simple model for the error in the measurement 
of the galaxy re dshift distributio n q(z)dz we used the convolution 
with a Gaussian JMa et al.ll2006l) . 



r(.z) 



1 



V2™1 



exp 



z 
'2al 



with (T- = 0.05, 



(49) 



which is of course too coarse to model errors in a proper weak lens- 
ing measurement, but will serve as an example. In contrast to the 
previous example, systematical errors in the parameters increase 
with larger numbers of polynomials used, and reach values of the 
order of the statistical error, in particular for n^ and erg, and to a 
lesser extend h. Quite similarly, the direction of the bias move the 
best fit point quickly away from the true cosmology because it is at 
roughly right angles relative to the degeneracy direction. 

Given the magnitude of systematical errors in comparison 
to their statistical accuracies, parameter estimation biases due to 
a wrongly chosen cosmology appear unlikely to affect measure- 
ments. In particular, if values estimated from unweighted tomogra- 
phy are used, the cosmological parameters are known well enough 
that TaRDiS-polynomials can be constructed in a reliable way and 
consequent parameter estimation biases would always be subdom- 
inant in relation to the statistical errors. In addition, it is always 
possible to iterate between parameter estimation and polynomial 
construction for narrowing down estimation biases. 



5.3 cross validation 

A possible way of validating the correctness of the assumed cos- 
mology for the construction of the orthogonal set of polynomi- 
als would be to estimate the cross spectra Cjj(i) which should 
be c ompatible with zero for the correct choice (similarly to 
iHute rer & White 2005), or alternatively, one can take advantage of 
the non-commutativity [5 , A'] = S N- NS 7^ of the signal covari- 
ance S and the noise covariance A'. If the polynomials have been 
constructed for the true cosmology with the first normalisation vari- 
ant eqn. l |23l l, 5 for unit-normalised polynomials is equal to the unit 
matrix and would always commute with N. In the case of the wrong 
cosmology, S and A' are symmetric matrices with different eigen- 
systems, and [5, Af] does not vanish. In Fig.[T4]this commutator is 
given as a function of multipole order C and depending on the num- 
ber of polynomials used. For giving a single number quantifying 
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Figure 14. Commutator ytr [S,N] between the signal covariance Sjj and 
the noise covariance Njj, cumulative in polynomial order q (indicated by 
colour), for the two systematics considered here: assuming a time-varying 
instead of a constant vvCDM model (wq = -0.8 and Wa = -0.2, solid lines) 
and choosing a too high matter density and a too high normalisation of P{k) 
(n,„ = 0.275 and erg = 0.85, dashed lines), both for EUCLID. 



the non-commutativity we chose to compute the trace tr [5, A']^, be- 
cause tr [S, A'] without squaring always vanishes. The stronger sys- 
tematic (detuning Q.„, and erg from their fiducial values) generates 
larger values compared to the wrong assumption of a time- varying 
equation of state of the dark energy component. Naturally, when 
using a single polynomial q = the commutator is always zero. 

Given the minor impact on choosing a wrong cosmologi- 
cal impact on parameter estimation from weak lensing spectra 
weighted with orthogonal polynomials, we point out that it should 
be possible to employ the polynomials iteratively by alternating be- 
tween parameter estimation and polynomial construction. At the 
same time we emphasise the possible usefulness of measuring cross 
spectra Cij{(), i t j, or the commutator tr [5, W]" for validating the 
cosmological model used for data analysis. 



6 SUMMARY 

The topic of this paper is a novel method of carrying out tomo- 
graphic weak lensing measurements by line of sight-weighting of 
the weak lensing convergence with specifically constructed orthog- 
onal polynomials. 

(i) The TaRDiS-polynomials have been constructed with a 
Gram-Schmidt orthogonalisation procedure in order to diagonalise 
the signal covariance matrix, and to project out statistically inde- 
pendent information of the weak lensing field by performing weak 
lensing tomography. They differ from traditional tomography and 
3-dimensional weak lensing methods in the respect that the signal 
covariance is diagonal instead of the noise covariance. 

(ii) The statistical error forcasts for cosmological parameters 
was investigated with the Fisher-matrix formalism. Because of the 
polynomial's property of providing statistically independent infor- 
mation of the weak convergence field, one sees a particularly sim- 
ple square-root scaling of the signal to noise-ratio and the statistical 
errors with the maximum multipole considered and the number of 
polynomials used. Clearly, statistical errors can be reduced when 
combining polynomials, where a small advantage can be observed 
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Figure 11. Parameter estimation biases (6^,5-,) in the parameters f2,„, o"8, h, n, and w, superimposed on the Icr-confidence regions, if the polynomials Piix) 
have been constructed for wCDM with an evolving equation of state parameterised by wq = -0.8 and Wo = -0.2 instead of wCDM with a constant equation of 
state with w = -0.9. In the estimation biases, the dot colour is proportional to the cumulative order q of the polynomials, and is plotted again as a function of 
q, for Hm (dots), erg (squares), h (lozenges), n, (triangles, pointing up), w (triangles, pointing down). For all computations, the EUCLID survey characteristics 



were used, and the data was accumulated up to f„, 



1000. 



in comparison to traditional tomograpiiy. Tiiis is due to tiie fact that 
in case of a matciiing between the assumed and the true cosmology, 
the Fisher-matrix assumes the largest possible values because of 
teh diagonalised covariance and provides consequently the small- 
est statistical errors. 

(iii) We extended the Fisher-matrix formalism for investigating 
how the assumption of a wrong cosmology in the construction of 
the polynomials impacts on the estimation of cosmological param- 
eters and to what extend biases are introduced. We assumed three 
systematically wrong priors: assumtion of a time varying instead 
of a constant dark energy equation of state parameter, as an exam- 
ple of a degree of freedom not contained in the model, a wrong £!,„- 
and CTg-pair, as being the two most prominent parameters determin- 
ing the strength of the weak lensing signal, and finally a convolved 
galaxy redshift distribution. All three systematics had a minor im- 
pact on the parameter estimation, with the biases decreasing if a 
larger number of polynomials was used. With 10 polynomials bi- 
ases were of the order of the statistical error even for strong mis- 
matches in the choice of the initial cosmology. Only the assumption 
of a wrong redshift distribution generates roughly constant biases, 
which, when normalised to the statistical errors, increase with the 
number of polynomials used. Our extension of the Fisher-matrix 



formalism treats the parameter estimation biases in full generality 
and does not assume a diagonal shape of neither the signal nor the 
noise covariance. 

(iv) The accuracy needed for constructing viable polynomials 
corresponds to the statistical error reachable with a simple un- 
weighted convergence spectrum, which is readily available, and can 
be improved by adding additional priors in the form of the CMB- 
likelihood, parameter constraints form baryon acoustic oscillations 
or from supernovae. 

(v) Additionally, it is possible to compute diagnostics for the 
choice of the prior cosmological model used for constructing the 
cosmology: Examples include the cross-spectra Sij(C) for two dif- 
ferent polynomials / t j, which should vanish for correctly- 
constructed polynomials, and the commutator [5, A'] = SN - NS 
between the signal and the noise covariance, which should vanish if 
the signal covariance is equal to the unit matrix for a certain choice 
of normalisation for the polynomials. 



We plan to extend our research on the usage of orthogonal 
polynomials to the weak lensing bispectrum in a future paper, and 
to carry out forecasts on dark energy cosmologies from combined 
constraints with spectrum and bispectrum tomography. 
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Figure 12. Parameter estimation biases (iS^.i?^) in tlie cosmological parameters n„,, erg, h, n, and w, superimposed on the Icr-confidence regions, if in the 
construction of the polynomials f2„, and erg were set too high (0.275 and 0.85 instead of 0.25 and 0.8, respectively) relative to the fiducial cosmology. The 
colour of the dots indicates the cumulative polynomial order q, for Q,,, (dots), erg (squares), h (lozenges), n, (triangles, pointing up) and w (triangles, pointing 
dowfn). For the plot, EUCLID survey characteristics were assumed, and the signal integrated up to (mux = 1000. 
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mials have been constructed for a linear (solid line) and a nonlinear (dashed 
line) CDM spectrum. 
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APPENDIX A: WEIGHTED LENSING SPECTRA 

Pig. lAll provides an alternative representation for Rig. |5] of the 
polynomial-weighted weak lensing spectra 5„(f). Here, the spec- 
tra are normalised to the value 5oo(^ = 10^)> such that the differ- 
ence is shape is easier to see. In particular the low multipoles are 
suppressed in amplitude by the weighting with Piix). 
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